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We study natural composite cold dark matter candidates which are pseudo Nambu-Goldstone 
bosons (pNGB) in models of dynamical electroweak symmetry breaking. Some of these can have 
a significant thermal relic abundance, while others must be mainly asymmetric dark matter. By 
considering the thermal abundance alone we find a lower bound of mw on the pNGB mass when 
the (composite) Higgs is heavier than 115 GeV. Being pNGBs, the dark matter candidates are in 
general light enough to be produced at the LHC. 



I. INTRODUCTION 



A second strongly coupled sector in Nature akin to 
QCD is a likely possibility. The new strong interaction 
may naturally break the electroweak (EW) symmetry 
through the formation of a chiral condensate, thus mak- 
ing the Standard Model (SM) Higgs a composite particle. 
Models of this type are called 'technicolor' (TC) [1 [1] 
and several new variants have been proposed recently 
[sHllI with interesting dynamics relevant for collid er p he- 
nomenology |12h15| as well as cosmology [13, fla - ISq . A 
review of these models and the phase diagram of strongly 
coupled theories can be found in Ref. [S^- A relevant 
point is that the technicolor dynamics is strongly modi- 
fied by the new interactions necessary to give masses to 
SM fermions [3g] and the interplay between these two sec- 
tors leads to an entirely new class of models, constraints 
on which were discussed in Ref. [39j . 

We discuss different possibilities for dark matter (DM) 
candidates within this rich framework and show that 
some of these composite states can be thermal relics while 
being sufficiently light to be produced at the LHC. 

We call dark matter candidates composed of tech- 
nicolor fields 'technicolor interacting massive particles' 
(TIMPs) and focus on those which are pseudo Nambu- 
Goldstone bosons (pNGB). Our analysis is general since 
we use a low energy effective description for the TIMPS 
which can easily be adapted for specific models. We will 
discuss some of these models [1, [ifl, ISfl ISOl which provide 
particularly interesting candidates for dark matter. 



II. THE SIMPLEST TIMPS FROM PARTIALLY 
GAUGED TECHNICOLOR 



An interesting class of TIMPs arise from partially 
gauged technicolor models [y, |40j in which only part of 
the TC group is gauged under the EW interactions. The 
EW gauged technifermions are organized in doublets in 
the usual way while the other technifermions are collec- 
tively denoted A-* , with / counting these flavors only: 
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These models were introduced originally in order to yield 
the smallest naive EW S parameter, while still being able 
to achieve walking dynamics. ^ The non-minimal flavor 
symmetry of the resulting model allows for a number of 
light states accessible at colliders. A similar scenario is 
envisaged in so-called conjormal technicolor J43l |4J] . The 
technifermions not gauged under the EW interactions es- 
sentially constitute a strongly interacting hidden sector. 

To be speciflc we consider a scalar TIMP, ~ AA, 
made of the SM gauge singlet technifermions A* and pos- 
sessing a global C/(l) symmetry protecting the lightest 
state against decay. Moreover we take (/> to be a pNGB 
from the breaking of chiral symmetry in the hidden sec- 
tor, which leaves this U{\) unbroken. This constitutes 
the simplest type of TIMP from the point of view of 
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^ The naive 5-parameter from a loop of technifermions counts the 
number of fermion doublets transforming under weak S'C/(2)l, 
while walking dynamics is required to reduce non-perturbative 
contributions to the full 5-parameter. The naive S-parameter 
has recently been conjectured to be the absolute lower bound of 
the full S-parameter |41|. 14211 . making the TC models presented 
here optimal with respect to satisfying the LEP precision data. 



its interactions so we study this first and consider later 
TIMPs with constituents charged under the SM. An ex- 
plicit model of partially gauged technicolor featuring this 
kind of TIMP is 'uhra minimal technicolor' (UMT) ^. 
We therefore identify our DM candidate with a complex 
scalar 0, singlet under SM interactions and charged un- 
der a new U{1) symmetry (not the usual technibaryon 
symmetry), which makes it stable. 

In addition to the TIMP we consider a light (com- 
posite) Higgs boson. A general effective Lagrangian to 
describe this situation is presented below and can be de- 
rived, for any specific model, from the UMT Lagrangian 
3 . At low energies we can describe the interactions of (j) 
through a chiral Lagrangian. The (composite) Higgs H 
couples to the TIMP as: 
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The interactions between technihadrons such as <j> made 
of EW singlet constituents and states with EW charged 
constituents (e.g. H) are due mainly to TC dynamics 
and, as such, the couplings between these two sectors are 
not suppressed [Sj. However, since (/) is a pNGB it must 
have either derivative couplings or the non-derivative 
couplings must vanish in the limit m^ — ^ 0. The mass 
m^ is assumed to come from interactions beyond the TC 
sector, e.g. 'extended technicolor' (ETC) [4^ |4g| which 
can provide masses for the TC Nambu-Goldstone bosons, 
as well as for SM fermions. The couplings di, ...,di are 
dimensionless and expected to be of 0(1), while A is the 
scale A ~ 47r_FV below which the derivative expansion is 
sensible. 

We emphasize the differences between a composite 
scalar TIMP and fundamental scalar dark matter con- 
sidered earlier (47-49]: i) The C/(l) is natural, i.e. it is 
identified with a global symmetry (not necessarily the 
technibaryon one), ii) Its pNGB nature makes the DM 
candidate naturally light with respect to the EW scale 
and influences the structure of its couplings, iii) Com- 
positeness requires the presence around the EW scale of 
spin-1 resonances in addition to the TIMP and the (com- 
posite) Higgs: their interplay can lead to striking collider 
signatures |10l.l30j. 



III. THERMAL VERSUS ASYMMETRIC DARK 
MATTER 

When the TIMP is a composite state made of particles 
charged under the EW interactions it becomes a good 
candidate for asymmetric dark matter, i.e. its present 
abundance is due to a relic asymmetry between the parti- 
cle and its antiparticle, just as for baryons. This has been 
the case usually considered when discussing TC DM can- 
didates [la, [S^ since the technibaryon self-annihilation 



cross-section, obtained by scaling the proton-antiproton 
annihilation cross-section up to the EW scale, is high 
enough to essentially erase any symmetric thermal relic 
abundance. Hence an asymmetry between technibaryons 
and anti-technibaryons is invoked, especially as this can 
be generated quite naturally in the same manner as 
for baryons. However, scaling up the proton-antiproton 
annihilation cross-section is not applicable to generic 
TIMPs, in particular not to pNGBs, hence they may have 
an interesting symmetric (thermal) relic abundance. 

Let us solve for the thermal relic abundance of TIMPs 
4> with singlet constituents using the Boltzman continuity 
equation |5ll.[52|: 
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where R is the cosmological scale-factor, n^ the TIMP 
number density and iTann the TIMP-antiTIMP annihila- 
tion cross-section (given in Appendix |X] along with the 
relevant interaction vertices). From the Lagrangian ^ 
we see that annihilations proceed via the (composite) 
Higgs into SM fermions and gauge bosons pairs, as well as 
into a pair of (composite) Higgs particles. As discussed 
in Refs. |53l - l55| . we can rewrite the continuity equation 
dS]) in terms of the dimensionless quantities Y = n^/s, 
Y^^ = n^^/s, and x = m^/T, where s = gsT^ is the spe- 
cific entropy determining the value of the adiabat RT: 
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and gp = p/T^ counts the number of relativistic degrees 
of freedom contributing to the energy density, which de- 
termines the Hubble expansion rate R/R. The values 
of gp{T) and gs{T) have been computed in the SM [5g| 
and are modified to account for the additional particle 
content of TC models. 

In the hot early universe, the particle abundance ini- 
tially tracks its equilibrium value but when the tempera- 
ture falls below its mass and it becomes non-relativistic, 
its equilibrium abundance falls exponentially due to the 
Boltzmann factor. Hence so does the annihilation rate, 
eventually becoming sufficiently small that the (comov- 
ing) particle abundance becomes constant. Defining the 
parameter A = (F — Y'^'^)/Y^'^, the freeze-out tempera- 
ture is given by |52]: 



Xfr^f ln[Aft.(2 + Af,)5,] 
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2b, 

1/3 



■In 



f ln(Af,(2-|-Afr)5,) 



(*)'. 



(27r)3/2"'' 



b-:"^\ 



(5) 



and g counts the internal degrees of freedom, e.g. g — 
2 for the TIMP. This gives a good match to the exact 
numerical solution of Eq.Q for the choice Af^ = 1.5 
which corresponds to the epoch when the annihilation 



rate, n!/'(craim v) , equals the logarithmic rate of change of 
the particle abundance itself: d In n'^^/dt = Xh R/R [5l| . 
Note that the usual criterion of equating the annihilation 
rate to the Hubble expansion rate R/R woul d gi ve an 
erroneous answer when there is an asymmetry [53| . 

We calculate the TIMP freeze out parameter ccfr as 
a function of the TIMP mass m^ taking A = 1 TeV, 
for three values of the (composite) Higgs mass tuh = 
250, 500, 1000 GeV. We also take the dimensionless effec- 
tive couplings to the (composite) Higgs to be of 0(1) and 
define di2 = di + ^2, ^34 = ^3 + ^4, since at low energies 
the di and ^2 terms contribute very nearly equally, as 
do the ^3 and d^ terms. There are spikes in x^ at the 
(composite) Higgs resonance when 2m^ = mu, however 
the simple a ppr oximation above is not reliable near such 
a resonance [57| and we must then solve the full continu- 
ity equation including the (composite) Higgs width. We 
do this using the programme MicrOMEGAs |58l460l| which 
computes the full annihilation cross-section of the model 
using CalcHEP [6l|. We also use LanHEP [12] for the model 
implementation. 

After freeze-out, only annihilations are important since 
the temperature is now too low for the inverse creations 
to proceed; the asymptotic abundance is then: 
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The resulting cosmological energy density of relic TIMPs 
is shown in Fig. [1] We have checked explicitly with the 
numerical code that the contributions from di and d2 
terms are (very nearly) identical, as are the contributions 
from d^ and (i4 terms. 
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FIG. 1: The relic TIMP {(f>) abundance vs. its mass. The 
thick lines show the Micromegas computation taking into ac- 
count the (composite) Higgs decay width. 

Fig. [5] shows the region in the (composite) Higgs ver- 
sus TIMP mass plane where the TIMP relic abundance 
matches the DM abundance fth'^ = 0.11 ± 0.01 (2cr) in- 
ferred from WMAP-7 [H. 

From Figs. [T] and [H we observe first that the relic energy 
density drops significantly for m^ ^ 2mH as expected 
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FIG. 2: Regions corresponding to Qh'^ = 0.11 ± 0.01 for the 
relic TIMP (0) abundance in the (composite) Higgs vs. TIMP 
mass plane. The dashed box shows that given tuh > 115 
GeV, we require m,j, > mw for TIMPs to be dark matter. 



(because (j) can then decay resonantly through the Higgs) 
and, second, that if mn is greater than about 115 GeV 
then for TIMPs to be dark matter requires m^ > mw- 

As discussed below, in the presence of an asymme- 
try the total relic abundance always increases relative to 
the same model with no asymmetry, so mw provides a 
general lower bound for the mass of the pNGB TIMPs 
we consider. It follows that in the interesting region 
mw ^ rrifi, < 1 TeV, symmetric relic TIMPs with singlet 
constituents could make up a significant fraction of the 
dark matter in the universe. However when the TIMP is 
heavier than about a TeV, the strength of the interaction 
is similar to that of an ordinary (scalar) technibaryon, 
and a relic abundance large enough to account for dark 
matter now does require an initial asymmetry similar to 
that of baryons as discussed earlier Id, 17]. We note that 
recently a different type of QCD-like pions (which do not 
carry a C/(I) quantum number) were also considered as 
dark matter candidates [63, KM . 



A. Adding an asymmetry 

To study the relic abundance in the presence of both a 
thermal component and an initial asymmetry we follow 
Ref.[53| and define the asymmetry as a = (1^+ — Y_)/2 
where Y± are the abundances of the majority and mi- 
nority species (TIMP and anti-TIMP) respectively. The 



abundance in thermal and chemical equilibrium is: 
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where fi is the chemical potential. The continuity equa- 
tion in the presence of an asymmetry is j53| : 
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\x-^ [y_-q(r^q + 2a) - r_(r_ + 2a)] , (8) 



and the total asymptotic abundance of TIMPs and anti- 
TIMPs is: 



n^h^ =5.5 X 10^(y_oo+a; 
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(9) 



In Fig. [3] we show the minority species abundance y_ as 
a function of x = m^/T for m^ = 100 GeV and a = 
9.8 X 10^^, 5.6 X 10~i°, 10"", taking the Higgs mass to 
be TTiH = 250, 500, 1000 GeV. 
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the relic abundance of (j) agrees with the DM abundance 
inferred from WMAP-7 [63| for a fixed value of the d 
coefficients. Just as in Figs. [1] and [5] we observe that if 
niH is greater than ~ 115 GeV we require m^ > niw to 
avoid an excessive TIMP relic abundance. However, due 
to the asymmetry, for TIMP masses above mw there is 
now a much broader range of mn and m^ where the relic 
abundance of d) matches the observed DM abundance. 
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FIG. 4: Regions in the (composite) Higgs versus TIMP mass 
plane corresponding to flh^ = 0.11 ± 0.01 for the TIMP {(jj) 
relic abundance, for different values of the relic asymmetry a. 
The dashed box shows that given niH > 115 GeV, we require 
m^ > mw for TIMPs to be dark matter. 



FIG. 3: TIMP ((/)) abundance when there is no asymmetry 
(thick lines), compared with the abundance of the minority 
species Y- (thin lines) when an asymmetry a is present. 



B. TIMPs with charged constituents 



We see from the figure that the symmetric component 
of the relic abundance is ~ 10% of the asymmetric com- 
ponent, when the asymmetry is comparable to the would- 
be symmetric abundance (in the absence of an asymme- 
try). In the limit where a <^ Y^o, the symmetric compo- 
nent is unchanged and the asymmetry provides a small 
addition to the total relic abundance. When a > Fqo, 
the abundance of the minority species is exponentially 
suppressed in a and provides a negligible addition to the 
asymmetric component. Adding an asymmetry will al- 
ways increase the relic abundance relative to its value in 
the same model in the absence of an asymmetry. This im- 
plies a non-trivial constraint on the mass of TIMPs with 
uncharged constituents discussed above, such as appear 
in e.g. the UMT model. The constraint is non-trivial 
since such states would evade direct detection at collid- 
ers unless the (composite) Higgs is very light [3^ , as well 
as the direct detection experiments discussed below. 

Fig. 21 shows the contour in the mn — m^ plane where 



We consider now pNGB TIMPs with charged con- 
stituents of the form T ~ UD arising from the TC sector 
carrying EW interactions (see Eq. [T]). These states carry 
an L'^(l) quantum number which makes them stable and it 
is natural to identify this global symmetry with the tech- 
nibaryon number. Such particles arise generally in TC 
models with the technifermions transforming in either 
real or pseudo-real representations of the gauge group. 
Explicit examples are furnished again in the UMT scheme 
in which the composite T is a SM singlet and the 'orthog- 
onal minimal technicolor' (OMT) model in which the T 
state is the isospin-0 component of a complex triplet [lOj . 

We demonstrate that similarly to the case of TIMPs 
with neutral constituents, TIMPs with charged con- 
stituents have a significant symmetric component in only 
a small region of parameter space. However, as opposed 
to the TIMPs with SM neutral constituents this region 
is essentially independent of the parameters of the (com- 
posite) Higgs interactions. In addition to the interactions 



investigated above, the scalar TIMPs containing charged 
constituents will also have an effective interaction with 
the photon, due to a non-zero electromagnetic charge ra- 
dius of T [il,[33l: 



lii 



ie-^T*blT d,F^' 



(10) 



The corresponding charge radius of the TIMP is tt ~ 
\fdBl t^- For our choice A = 1 TeV we consider the range 
Ids I = to 0.3 while for a higher cut-off A, a larger 
d-B ^ 0(1) would be expected. 

In case of the TIMP T there are also contact inter- 
actions with two SM vector bosons V ^ arising from the 
kinetic term of the chiral Lagrangian, which can signifi- 
cantly affect the symmetric relic density. In general these 
can be written as 
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vv = \t*TV^V^' Tx \ [As,[As,Xt]]Xt. (11) 



[As,[As,Xt]]Xt. 



where X^ is the generator of the broken (techni) flavor 
direction corresponding to the TIMP and As are the, 
appropriately normalized, EW generators imbedded in 
the TC chiral group ^^Q^. The resuhing TTWW and 
TTZZ contact interactions (assuming that the EW sym- 
metry is broken already) are: 
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■TT[dwW^W^' + dzZf,Z''] , (12) 



with dw = 5^ and dz = (ff^ + g''^)/2 for the TIMPs T of 
both the UMT and the OMT models. 

In Fig. [5] we display the effect of the W and Z con- 
tact interactions, as well as of the charge radius interac- 
tion, on the thermal relic abundance. It is seen that for 
rriT significantly below mw, the interactions in Eg. ljlOp 
do affect TIMP annihilations significantly, even so the 
TIMP relic abundance is too large. As rriT increases to- 
wards mw, annihilations due to the TTVV interactions 
in Eq. (fT2)) begin to dominate and reduce the TIMP relic 
abundance to the observationally acceptable level only 
when the TIMP is a few GeV lighter than the W. 

Again, once we include an asymmetry, the range of 
TIMP masses where a cosmologically acceptable sym- 
metric relic abundance is achieved, can be much broader. 
This is demonstrated in Fig. [S] which shows the contours 
of flh^ for different values of the asymmetry and rriT. 
There is a cross-over in the small TIMP mass region men- 
tioned above, from just below mw where the relic abun- 
dance is dominated by the asymmetry, to just above mw 
where it is dominated by the symmetric component. 
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FIG. 5: The TIMP (T) relic energy density as a function 
of its mass when W and Z contact interactions are present, 
calculated taking into account their ofF-shell decays. The hor- 
izontal band corresponds to Qh^ = 0.11 ± 0.01. 
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FIG. 6: The region in the asymmetry (a) vs. TIMP (T) mass 
plane which yields the marked relic energy density (consistent 
with WMAP) when W and Z contact interactions are present. 



IV. DIRECT DETECTION 

Direct detection of the TIMPs (j) with SM singlet con- 
stituents will be challenging. The exchange of the (com- 



posite) Higgs leading to a scattering cross-section on nu- 
clei is the most relevant interaction here, and we will as- 
sume that the (composite) Higgs couples to SM fermions 
with ordinary Yukawa couplings. (In fact this represents 
an upper bound and so the actual cross-section could be 



lower). Again, since the TIMP is a pNGB the couphngs 
to the Higgs are suppressed at low masses. The TIMP nu- 
cleon scattering cross-section from the (composite) Higgs 
exchange is given by: 
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dH = {di+d2)j^, (13) 



where fi is the nucleon-TIMP reduced mass, v the elec- 
troweak vev and / parameterizes the (composite) Higgs- 
nucleon coupling. We refer to Refs.|73. uM ^^^ recent 
discussions on the strange quark contribution to / which 
we take to be / = 0.3. 

For TIMPs T with charged constituents there is an 
additional contribution to the scattering on nuclei via 
the charge radius operator [la . |30| 
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To take into account the possible interference between 
the composite Higgs and photon exchange [30| we write 
an averaged scalar cross-section per nucleon as 
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The direct detection cross-section per nucleon as a func- 
tion of the TIMP mass is shown in Fig. [3 where we also 
indicate, (following Ref . (7411 ) the limits from the CDMS 
II [zi] and XENON-100 [731 experiments. 

For Higgs exchange only {ds = 0) the direct detec- 
tion cross-section increases with TIMP mass since the 
TIMPs are pNGBs (c./. Ref. iS^J where the TIMPs were 
not derivatively coupled to the (composite) Higgs). How- 
ever, the charge radius interaction can significantly al- 
ter the cross-section, by up to 2 orders of magnitude 
for I^bI ~ 0.3, thus greatly affecting the discovery po- 
tential of direct detection experiments. For example, for 
de = -0.3 and di2 = 3 the signals from the TIMP T with 
SM charged constituents would have been observed al- 
ready by CDMS II and XENON-100. On the other hand, 
for positive values of ds, there is a destructive interfer- 
ence between (composite) Higgs and photon exchange, 
lowering the cross section to 10^*^ cm^ at tjit — 200 GeV 
for a particular choice of the parameters. 

Since the direct detection rate depends strongly on the 
(composite) Higgs mass, we present in Fig. [5] results in 
the {Mh—ttit) plane for different values of the ds param- 
eter. The figure also shows the exclusion limits from the 
XENON-10 [zi], CDMS II 76] and XENON-100 [T^ ex- 
periments. One can see that for d^ = and di2 = 3 (top 
frame), XENON-100 and CDMS II can cover essentially 
the whole range of rriT for Mh below 150 GeV where the 
cross-section always exceeds 10~^^ cm^. Our results triv- 
ially scale as d^^ for ds = 0. Negative ds = —0.3 signif- 
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FIG. 7: Direct detection cross-section (per nucleon) for TIMP 
scattering off nuclei. The full line is for Higgs exchange only 
(with rriH ~ 200 GeV) while the short- and long-dashed lines 
show the additional effect of the charge radius operator (with 
cIb = —0.3 and dg = -fO.3 respectively). The shaded region 
is experimentally excluded. 



icantly enhances the cross-section (middle frame) while 
positive ds = 0.3 (bottom frame) brings the negative 
interference effect into play, resulting in a deep valley. 



V. CONCLUSIONS 

We have shown that models of dynamical EW sym- 
metry breaking can provide symmetric (thermal) DM, as 
well as asymmetric (non-thermal) DM. This is true in 
particular for partially gauged technicolor [y, |40[ which 
can satisfy constraints from EW precision measurements. 

From our analysis we conclude that for pNGB TIMPs: 

1) The TIMP cannot be significantly lighter than the 
W if its relic abundance is to be acceptable, unless the 
(composite) Higgs mass is below 115 GeV — this holds 
whether it is symmetric, asymmetric or a combination. 

2) If the TIMP is made of constituents charged under 
EW interactions, it can be symmetric DM only in a nar- 
row mass range close to mw- Above this mass an initial 
asymmetry is required for TIMPs to be dark matter. 

3) If the constituents of the TIMP are neutral with 
respect to the EW interactions then it can be symmetric 
dark matter for a range of masses tied to the (composite) 
Higgs mass. This does not exclude the possibility that it 
has an asymmetry as well. 

4) Direct detection of light pNGB TIMPs with neu- 
tral constituents is challenging due to its mass sup- 
pressed couplings to the (composite) Higgs. However, 
for TIMPs with charged constituents there is an addi- 
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FIG. 8: Contours of the direct detection cross-section (per 
nucleon) for TIMP scattering off nuclei in the TIMP mass vs. 
(composite) Higgs mass plane. The top, middle and bottom 
frames are for ds = 0, —0.3, +0.3 respectively. The dashed 
curves show the upper limits from CDMS II (red), XENON- 
10 (blue) and XENON-100 (green). 



tional charge radius interaction which, if sizeable, can 
bring such TIMPs within the reach of current nuclear 
recoil detection experiments. 

Both types of TIMPs considered here - with charged 
or neutral constituents - can co-exist and contribute to 
the dark matter (as in e.g. the UMT model Q). These 
models provide interesting signals for direct dark matter 
detection experiments which are already sensitive enough 
to exclude TIMPs in certain regions of parameter space 
or even discover them in the near future. 
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Appendix A: Annihilation cross-section 

The relevant vertex factors at low energies (if we keep 
only the light (composite) Higgs in the spectrum) are: 



i>'(l>H 



*<j)HH 



A 



P4>*P4> 
A2 



di+ d-z- 



da + d4' 



A 



A2 



idi 



idsi 



A 



A2 



(Al) 



where di2 = rfi + ^2 and ^34 = ^3 + 0^4 are the only inde- 
pendent parameters at low energies. Hence the (compos- 
ite) Higgs mediated contributions to the cms annihilation 
cross-section (crwroi) in the limit Wrci — >■ are: 
dxh* -^ HH : 
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Here the fermion Yukawa coupling is A/ = rnj /v where 
V ~ 246 GeV and mf is the fermion mass, while c/ = 1,3 
for leptons and quarks respectively. The contributions 
from the photon mediated annihilations from the charge- 
radius operator are negligible and we do not include 
these. 



We implement the Lagrangian in Eq. ^ in CalcHEP 
eil and in MicrOMEGAs [5^ (using the LanHEP module 



62| to check the above implementation) , in order to com- 
pute the full 2 — >■ 2 annihilation cross-section including 
finite widths and to study the collider phenomenology. 
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